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We present a complete calculation of the cross section for neutralino annihilation into 
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the two-gluon final state. This channel can be quite important for the phenomenology of 
neutralino annihilation due to the well-known helicity suppression of neutralino annihila- 
tion into light quarks and leptons. In addition, we calculate the cross section for annihi- 
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lation of neutralinos into a gluon and quark-antiquark pair, and discuss QCD corrections 
to the tree-level cross sections for neutralino annihilation into quarks. If the neutralino 
^ . is lighter than the top quark, the effect of these results on high-energy neutrino signals 
from neutralino annihilation in the Sun and in the Earth can be significant, especially 
if the neutralino is primarily gaugino. On the other hand, our results should have little 
effect on calculations of the cosmological abundance of neutralinos. We also briefly discuss 
implications for cosmic-ray antiprotons from neutralino annihilation in the galactic halo. 
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I. Introduction 

Luminous matter almost certainly does not account for all matter in the Universe [lj , 
and this matter deficit inspires both particle-physics and astrophysics speculation. Among 
the particle-physics solutions discussed, perhaps the most attractive idea is that stable 
weakly interacting massive particles (WIMPs) may make up the dark matter. Currently, 
the most promising candidate WIMP is the neutralino 0, a linear combination of the 
supersymmetric partners of the photon, Z , and Higgs bosons. 

It is well known that the cosmological abundance of a WIMP is inversely proportional 
to its annihilation cross section. In the earliest calculations, annihilation of neutralinos 
into light fermions was considered |§f§]; subsequently, annihilation into gauge and Higgs 
bosons was taken into account At this point, the cross sections for annihilation 

into all two-body final states which occur at tree level had been calculated. The basic 
conclusion of all these papers is that, in almost all regions of supersymmetric parameter 
space, the neutralino makes an excellent candidate for the dark matter in galactic halos, 
and, in some cases, can provide an abundance suitable to account for a flat Universe. 

The neutralino annihilation cross section is also needed to determine event rates in 
numerous schemes for indirect detection of WIMPs in the galactic halo. The existence of 
WIMPs in the halo may be inferred through observation of distinctive spectra of cosmic- 
ray antiprotons, 7-rays, or positrons produced by annihilation of WIMPs in the halo || . 
Perhaps a more promising method of detection involves observation of energetic- neutrinos 
from WIMP annihilation in the Sun and/or in the Earth P| ||10|| . If neutralinos reside in the 
galactic halo, then they will be captured in the Sun and in the Earth [11 1. Annihilation of 
neutralinos therein would produce, amongst other things, neutrinos whose energy is some 
fraction of the neutralino mass, which could be detected in current (e.g., 1MB, Kamiokande 

II, MACRO) or next-generation (e.g. DUMAND, AMANDA, super-Kamiokande) experi- 
ments fr^| . 

In this paper, we calculate the cross section for annihilation of neutralinos into two 
gluons (xx —* 99) and discuss the implications for indirect-detection searches, especially 
those involving observation of energetic-neutrinos from neutralino annihilation in the Sun 
and in the Earth. We also calculate the cross section for annihilation into the quark- 
antiquark-gluon final state (xX ~ * QQ9)i an d we improve the tree-level cross sections for 
annihilation of neutralinos into quarks by including leading-log QCD corrections. 

There are a few simple reasons to believe that the two-gluon final state may be im- 
portant for indirect-detection experiments, although we do not expect this final state to 
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have much impact on cosmological relic-abundance calculations. Neutralinos in the halo, 
Sun, and Earth move with velocities v < 10 -3 , so annihilation through P waves (i.e., an- 
gular momentum / = 1) is suppressed by roughly v 2 < 10 -6 ; in other words, annihilation 
can occur only through an S wave. Neutralinos are majorana particles, so the S-wave 



hand, there is no such suppression of S'-wave annihilation into gluons. Therefore, even 
though annihilation into gluons is formally suppressed relative to that into quarks by a 2 , 
the square of the strong coupling constant, in practice, S'-wave annihilation of neutralinos 
into gluons may be comparable to, or even stronger than annihilation into quarks, as first 
noted by Rudaz and Bergstrom ||. For given SUSY parameters, if the neutralino annihi- 
lates primarily into light fermions at tree-level, we expect our results to be important |ITJ . 
This will be the case if the neutralino is lighter than the W boson. It may also be the 
case if the neutralino is primarily gaugino and heavier than the W boson but lighter than 
the top quark || . In addition, since S-wave annihilation into light fermions is additionally 
suppressed in the limit of large neutralino mass the gg annihilation channel should become 
increasingly important as the neutralino mass is increased. 

When neutralinos freeze out in the early Universe, they are moving with velocities of 
order 0.5, so the cross section for annihilation contains significant P-wave as well as S-wave 
contributions. Since P-wave annihilation into light fermions is not suppressed in this case, 
annihilation into gluons is suppressed by a 2 relative to that into light fermions. Therefore, 
inclusion of the gg final state should have no more than a small effect on the cosmological 
abundance, and thus we do not consider these abundance calculations further. 

In the following Section, we describe the calculation and present results of the cross 
sections for XX 99 an d XX ~^ QQ9- We also discuss the leading-log QCD correction 
to the annihilation cross section for XX ~ * QQ- We present our results for a neutralino 
of arbitrary mass and mixing and for arbitrary squark masses and mixings. From these, 
the cross sections for any given supersymmetric parameters can be obtained. In Section 
3, we illustrate our results for a few simple examples and discuss the importance of the 
new channels in various regions of parameter space. In Section 4, we review production 
of energetic neutrinos from annihilation of neutralinos in the Earth, and we discuss the 
effect of the new cross sections on the predicted event rates. In the final Section, we 
briefly summarize and comment on the implications of our results for predicted fluxes of 
cosmic-ray antiprotons produced by annihilation of neutralinos in the halo. 
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2. Annihilation Cross Sections 

In this Section, we review some relevant supersymmetry (SUSY) formalism and 
present results for the cross sections. We use the conventions and notation of Ref. 0. 
There are four neutralinos which are linear combinations of B, W3, h\, and h®, the super- 
symmetric partners of the U(l) gauge field, the third component of the 577(2) gauge field, 
and Higgs fields, respectively In the (S, TU3, h®, h®) basis, the neutralino mass matrix is 

/ Mi — mzsw cos (3 mzswsm(3 \ 

M 2 mzcw cos (3 —mzcwsin(3 , . 

—mzsw cos f3 mzcw cos (3 — fi ' 

V mzSw^inP —mzcw sin (3 — [i / 

where \i is the higgsino mass parameter, M\ and Mi are the gaugino mass parameters, 
sw = sin6*iy, and cw = cos$w, and we adopt the GUT relation, M 2 = (5/3)Mi tan 2 9w 
The ratio of Higgs vacuum expectation values is tan (3 = {H^) / The elements of 

the matrix that diagonalizes the mass matrix are iVjj. Unlike in Ref. 0, we allow the 
neutralino mass eigenvalues to take on negative values, so the Nij are always real. We 
denote the lightest neutralino by the subscript "0" ; specifically, the lightest neutralino is 

X = N 01 B + N02W3 + N03M + N 4~h 2 . (2.2) 

2.1. Cross section for XX ~^ 99 

In this subsection we present the results of our calculation of the cross-section for 
annihilation of non-relativistic neutralinos into two gluons. Our results are only for the 
limit of vanishing relative incident velocity, v — > 0, and are therefore suitable for use in 
indirect-detection calculations. Although they are not valid for relic-abundance calcula- 
tions, they could be used to estimate the magnitude of the annihilation cross section at 
v ~ 0.5. 

In Fig. 1, the one-loop diagrams for annihilation of neutralinos into two gluons are 
displayed, and anti-symmetrization (symmetrization) on the identical particles in the initial 
(final) state is understood. Due to the helicity structure in the v — > limit, only diagrams 
(a), (b), (e), and (f) contribute. This is because the CP eigenvalue of the initial two- 
neutralino state in this limit is —1. The two gluons in the final state could therefore be 
produced only by an operator of the form G^G^ , where G a is the gluon field-strength 
tensor. Inspection of the Dirac structure of diagrams (c) and (d) shows that they cannot 
give rise to such a term. Furthermore, the Dirac structure of the remaining diagrams also 
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e. f. 

Fig. 1 Feynman diagrams for annihilation of neutralinos into two gluons. Figures (a)-(d) 
show quark-squark loops, while (e) and (f) show the exchange of the Z° boson or the 
pseudoscalar Higgs boson A . 

simplifies considerably in the v — > limit. This allows us to write a relatively compact 
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form for the final result. 

The matrix element for XX ~~ ^ 99 is 



M = -^e(k ll k 2 ,e u e 2 )(-l) h+1 H hrh M^8 ab . 



(2.3) 



The imaginary part of M. is 
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where /3 9 = y^l — m 2 /m x 2 . The real part of M. is 



_ / [ — ln 

2m x Jo V x 



x 2 a + x(b - 1 - a) + 1 



+ 
+ 



S g b + D a Vabf I 1 



1 + a - 6 
1 



1 — x 1 + x 



1 — 6 + xa 
x ln 



S„b ( - + 

x 



1 



1 -x 



ln 



+ D a \fab 



-x 2 a + x(b- 1 + a) + 1 



x 2 a + x(6 - a - 1) + 1 



6 + a(l - x 2 ) 
1 



1 -x 



x 2 a — x(a + b — 1) — 1 



+ 



1 



6 — 1 + ax 
x ln 



S q b 



1 



1 



1 + x 



- D q Vab — — - 

X + 1 



x 2 a + x(b — 1 — a) + 1 

2 



+ 21 



m q ( 2 9A qq 9Axx 



+ 



x 2 J m x \ (4m x 2 — m\) ' m x 3<3 m;| cos 6>jy 



9zxx9 



(2.5) 



Note that there is a sum over quarks, and for each quark, there is an additional sum 
over the two squarks in the term containing the parametric integral. We have defined 
S q = A 2 + B 2 , and D q = A 2 — B 2 , and A q and B q are quark-squark-neutralino couplings 
defined below. We have also defined a = m x 2 /m~ and b = m 2 /m~, and it is important to 



note that the sign of v a& is the sign of 



Generally, left and right squarks (q~L and q R ) may mix [H], and the squark mass 
eigenstates are then (see, e.g., Ref. 0), 

Qi = q~L cos 9 q + q R sind q , 

(2.6) 

<?2 = -Ql sm 6> ? + q R cos 6> 9 , 

where 6 q are the squark mixing angles. The squark-quark-neutralino couplings for the 
lighter squark eigenstates, denoted by the subscript "1" , are 

A q = -[cosflgpfg + Z q ) + sm6 q (Y q + Z q )}, 

I (2-7) 
B q = -[cos6 q {X q - Z q ) + sin6 q (Z q - Y q )). 

The squark-quark-neutralino couplings for the heavier squark eigenstate, denoted by a 
subscript "2", are obtained by making the substitutions sin6> g — > cos6> g and cos^ — > 
— sin^ g . Here, we have defined 

X q = -V2[gT 3q N 2 - g'(T 3q - e q )N 01 ], 

r , (2 - 8) 
Y q = V2g'N 01 e q . 

For up-type quarks we have 

\J imw sin p 

and for down-type quarks we have 

Z q = 1. q ° — . (2.10) 

V2mpycos,9 

In the above, T 3q = ±1/2 is the weak isospin of the quark, e q is its electric charge, g is the 
SU(2)l gauge coupling, g' is the U(1)y gauge coupling, and mw is the VF-boson mass. 

The mass of the pseudoscalar Higgs boson, A (sometimes referred to in the literature 
as H®), is ttia; the couplings of the pseudoscalar Higgs boson to up- type quarks are 

gm q cot (3 

9A «« = (2 - n) 

and the couplings to down-type quarks are obtained by making the substitution cot (3 — > 
tan j3. The coupling of the pseudoscalar Higgs boson to the neutralinos is 



9Axx = -(N 02 -tcm9 w N 01 )(N 03 smf3-N 04 cosP), 
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(2.12) 



and the coupling of the Z° boson to the neutralinos is 



9Zxx 



t2 
'03 



4 cos 9w 



(2.13) 



The function I(x), which arises from the three-point function in diagrams (e) and (f) 
of Fig. 1, is given by 



I(x) = 



i [>- 2 m - - 2 

l 2 

arctan J— = 



for x < 1, 
for x > 1, 



(2.14) 



where /3(x) = Vl — x - 

Note that the integration over the Feynman parameter x can be transformed in many 
ways. In particular, by judicious use of the transformation x —>■ 1 — x in some of the terms 
of the integrand, it is possible to improve the convergence of numerical integration. A 
poor transformation of this type would yield singularities at the endpoints which would, 
in principle, cancel when integrated. As we have written it above, the integrand has been 
transformed so that it has no singularities at the endpoints, and this is the form most 
suitable for numerical integration. 

Given the matrix element above, the cross section times relative velocity is 



°99 V 



■ a s m X 

8tt 3 



(2.15) 



Although the calculation is lengthy and complicated, our result for the cross section 
was obtained independently by several of the authors. In addition, the diagrams for neu- 
tralino annihilation into two gluons are similar to some that appear in the calculation for 
annihilation into two photons. Rudaz performed the calculation for pure photinos and 
higgsinos in the limit of large squark masses, Giudice and Griest obtained similar results 
and generalized to neutralinos of arbitrary mixing, and Bergstrom performed the calcu- 
lation for photinos for arbitrary squark and photino masses ||. For large squark masses 
(a, 6< 1), the expression for the real part of the matrix element, Eq. (|2.5|) , simplifies to 



ReM 




mi ^-^ 



+ 21 



m 



i ■! 
2 



-■% + ( 



m 



x 



m<> 



m r 



m< 



m. 



J2 

X 



2Wg 9Aqq9AxX 

m x (4m x 2 — m\) 



+ 



m r 



m x 2 m 2 z 



hq9zxx 



9 



cos 9 



w 



(2.16) 




Fig. 2 Feynman diagrams for the tree-level annihilation of neutralinos into a gluon and a 
quark antiquark pair. 



If the neutralino is a pure photino, the terms due to Z° and A exchange disappear, and 
D q —* 0. In this limit we reproduce Bergstrom's result || for the box diagram numerically, 
but our expression (|2.5| ) is significantly more compact, involving only a single (rather than 
double) integral over Feynman parameters. Note also that, even in the limit m q — > 0, 
the cross section is nonzero which reflects the fact that there is no chirality suppression of 
<S-wave neutralino annihilation into gluons. 

2.2. Cross section for XX ~^ 9QQ 

In this subsection we give the cross-section for the annihilation of neutralinos into a 
gluon and quark-antiquark pair. This cross section also remains finite as m q — > 0, as does 
the gluon-pair cross section discussed in the previous subsection. Moreover, it contains 
only a single strong-interaction vertex, i.e. only one factor of a s . If the neutralino is 
heavier than the top quark, annihilation into top quarks is unsuppressed, and the gqq 
final state will be negligible in comparison. Therefore, we consider only the case that the 
neutralino is lighter than the top quark, and calculate this cross section in the limit of zero 
mass for the final state quarks. 

In Fig. 2, the diagrams for the cross section for XX ~^ 9QQ are displayed. It is worth 
noting that the expected infra-red divergence from soft-gluon emission in these diagrams 
does not occur in the zero quark-mass limit. Diagrams (a) and (b) do not vanish in this 
limit, but their contributions are infra-red finite and combine with the infra-red finite 
contribution of diagram (c). As in Fig. 1, the anti-symmetrization on the identical initial 
particles in Fig. 2 is understood. Note that the contributions from Z° or A exchange in 
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the s-channel vanish in the massless-quark limit, so that the diagrams of Fig. 2 are the 
only ones to consider. 

After some algebra, we reduce the calculation to the following phase space integral; 



^ 8A*a a f 1 f 1 
> n o / dxi / dx 2 &(xi + x 2 - 1) 
V *" m x Jo Jo 



(2.17) 

(xi + x 2 - l)(xf + x%- 2(xi + x 2 - 1)) 
(l-2x 1 -r 9 -) 2 (l-2x 2 -r 9 -) 2 ' 

where Xi = Ei/m^ are the normalized energies of the quark and antiquark, = m~/m^ 2 , 
and A g is the quark-squark-neutralino coupling defined above (A q = ±B q for massless 
quarks). The displayed form of the integrand is most convenient as a representation for 
the double differential cross-section, do j dx\dx 2 . One of the remaining integrals in ( |2.17| ) 
is elementary; integrating this gives 



gqq Y 27r2m x 2 ^o (l + r 9 -)(2x + r g --l 

where 



fi(x) = —lAx — 12rqX — 6r~x + 18x 2 + 2rqX 2 — 8x' : 



q u, q ^ i j-wo, i ~, q 

"\ + 3r~ — 12a: — 16rqx ,., (/ 
The remaining integral can be treated numerically in direct fashion. 



f 2 {x) = 7 + 13rq + 9r? + 3r? - 12x - 16r g -x - 4r~x + 4x 2 + 4rqX 2 . 



(2.19) 



2.3. Cross Section for XX ff 

For clarity and comparison, we present the cross section for annihilation of neutralinos 
into fermions and quarks 00- In the limit of zero relative velocity, the cross section is 



c f 9(m^ 2 - m 2 ) 2 / m 2 

^ — ^ — V 1 ^7 

1 



x 



m-~ + m^ 2 — m 2 



(2.20) 



where the sum on / is over leptons as well as quarks, and the color factor Cf is 3 for quarks 
and 1 for leptons. This cross section is proportional to the square of the quark mass. 
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(Recall that \Af\ = \Bf\ for massless fermions, i.e. Df oc m/.) Also note that the cross 
section is proportional to the square of the imaginary part of the amplitude for neutralino 
annihilation into gluons, as it should be. 

Previously, the cross section for neutralino annihilation into quarks has been evaluated 
only at tree-level. In addition to the tree-level contributions, there will be QCD corrections 
to this cross section which can be important. These QCD corrections can broadly be 
grouped into three classes: Emission of a hard gluon at large angle to both quark and 
anti-quark; emission of a hard gluon (almost) collinear to either quark or anti-quark; and 
virtual corrections plus soft-gluon emission (these two must be added to cancel infrared 
divergences) . The first correction has been treated in the previous subsection, where we saw 
that it remains finite as m q — > 0. The main effect of the second contribution is to change 
the energy spectrum of the produced quarks; this will be included using Monte Carlo 
results, as explained in Sec. 4. Finally, the third class of corrections only renormalizes the 
overall value of the cross section. As usual, the leading logarithms in this correction can 
be most easily taken into account by introducing "running parameters" . This procedure, 
which is based on the renormalization group, also automatically re-sums all powers of these 
logarithms as they arise in higher orders of perturbation theory. 

In the case at hand there are two parameters that "run" due to QCD corrections: The 
Yukawa contribution Z q to the LSP-quark-squark couplings, and the quark mass m q . Eqs. 
(2.9), (2.10) show that Yukawa couplings and masses are proportional, and thus run in the 
same way. Running quark masses have previously been utilized in calculations of Higgs 
boson decay widths fl5|] , but they are also useful for the treatment of QCD corrections 
to the production of massive quark pairs in e + e~ annihilation [IB], which more closely 
resembles our case of LSP annihilation. 

For momentum transfer m b < Q < mt, where mj and m t are the bottom- and top- 
quark masses, the running masses are given by |15[ 



, , , . 12/23 

a, " 1 ■ ' 



m b (Q) = m h {m b ) ( ^9)\ ( 2 .21) 
\a s {m b ) J 



and 

,{mb)y" m ( <*s{Q) 



, , 12/25 / x 12/23 

CK C " /,></), \ 



m c (Q)=m c (m c ) "T^T 5 ( 2 - 22 ) 

a s {m c ) J \a s {m b ) 



if Q > m t , the last factor becomes 

frr > \ 12/23 / / n, x 12/23 / / \ 12/21 

a s (m b ) ) \a s (m b ) J \a s (m t ) 
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(2.23) 



where mf,(m&) = 4.5 GeV and m c (m c ) = 1.35 GeV are the on-shell quark masses, and 
a s (Q) is the strong coupling constant at scale Q [jIT| . We use Q = in our calculation. 

For neutralino masses in the range of O(10 — 100) GeV, the running of the quark 
masses is significant, and since the tree-level cross section is proportional to the square 
of the quark mass, this affects the cross sections for annihilation into quarks significantly. 
For example, if the neutralino mass is 80 GeV, the running c-quark mass is about 3/5 its 
tree- level value, and the running 6-quark mass is about 3/4 its tree- level value, so the cross 
section for annihilation into quarks is roughly half that suggested by using the tree-level 
masses. The consequences for the rate of energetic-neutrino production from neutralino 
annihilation will be discussed below. The running of the quark masses should have little 
effect on cosmological relic-neutralino abundances because the cross section for neutralino 
annihilation through a P wave is not governed by the quark mass. 

3. Illustrative Examples 

The cross sections for the various annihilation channels are complicated and depend on 
many parameters making it difficult to analytically assess the relative importance of each 
channel for arbitrary input parameters. Therefore, to illustrate the possible importance of 
the new channels, at least for a neutralino of some given simple composition, we consider 
several specific examples. 

The first limit we consider is that where the neutralino is a pure photino: iVoi = 
cos6>vk, N Q2 = sin 9w, and N 03 = N 04: = 0. In addition, we also assume that the light- 
fermion masses are small (m q <C "^x), an d that the photino mass is negligible compared 
with the squark and top-quark masses (m x <C m^mt). Furthermore, the pure photino 
has no coupling to the Z° or A bosons. We also take degenerate squarks and assume 
no squark mixing. Then, D q — and S q = 4ivae q , where a is the electromagnetic fine- 
structure constant. We also note that I(x) ~ 1/x for x 1. With these assumptions, the 
cross section for photino annihilation into gluons is 

2 

a gg v ~4a 2 o^^. (3.1) 

With the same assumptions, annihilation into fermions occurs primarily into fr pairs and, 
to a smaller degree, into cc, and bb pairs. The cross section for annihilation into these 
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channels is 

^-^ 8ircfa 2 eim 2 f 

a ff V ~ 1^ ^4 

/ /' , s 

, x 2 (3-2) 

120 a 2 ( 1 GeV 



We have not used running quark masses here. Use of running quark masses decreases 
this estimate only slightly due to the slow running of the r-lepton mass. With similar 
assumptions, the cross section for annihilation into the 3-body final states are suppressed 
by larger inverse powers of the squark mass. 

Comparing the annihilation cross sections for the pure photino, we find the following 
ratio: 

valid for nib *C <C nit, and <C rriq. Therefore, the relative importance of the gluon 
annihilation channel increases roughly with the square of the photino mass and becomes 
dominant for photino masses greater than about 45 GeV. 

We now consider the pure S-ino limit: iVoi = 1 and iVo2 = Nq% = N 04: = 0. If we use 
the GUT relation Mi = (5/3)M 2 tan 2 6 W and if M 2 <C n, then the neutralino is usually 
primarily S-ino. In fact, in most cases where the neutralino is primarily gaugino, it more 
closely resembles a S-ino than a photino. Once again, we take mj <C <C m t ,mq. In 
this case, the cross section for S-ino annihilation into gluons is 

49\ 2 g' 4 a 2 m- 2 



while the cross section for annihilation into light fermions is 

°n« = £ + Sf ? = 0.3 ^ (i^) 2 . (3.5, 



Again, using bare quark masses, we find the ratio 



^99_^ ( m x 



f f V 50 GeV 



(3.6) 



a result similar to that in the photino case. 

The next limit we consider is the pure-higgsino limit: iVoi = N 02 = 0, and N s = 
iVo4 = l/y/2. We make the same assumptions about the masses. Again, in this limit, 
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D q = 0. The main contribution to the gg cross section comes from t-quark loops, unless 
tan j3 is very large, and the cross section is roughly 

°99 V ~ 1B 3 4 V 3 - 7 

where we have assumed m^ 2 <C m 2 . For the light-fermion final states, annihilation occurs 
primarily into b quarks, and 



3S 2 m 2 h 
2nm 4 - 



'ff ^ "dbr- ( 3 - 8 ) 



Using on-shell quark masses we find that 



°99 _ ( m X 



) 6 cot 4 /3. (3.9) 



a ff V 20 GeV/ 

If we use running quark masses, the number is closer to 15 GeV. We caution, however, 
that the LSP is only very rarely sufficiently pure higgsino for Eq. ( [OP to be applicable 
even approximately. 

In the previous examples, we have considered neutralinos with no coupling to the 
Z° or A bosons. Now let us consider what happens if the squarks are heavy enough 
that annihilation into gluons and fermions occurs primarily through, for example, the A°- 
boson resonance. Since two neutralinos in an S'-wave cannot produce a physical Z° boson, 
the contribution from Z° exchange is almost always subdominant. Again, we consider 
nib rat- Then, decay of the intermediate virtual A occurs primarily through 

the top-quark loop. Of the final-state light fermions, decay occurs primarily into bb pairs. 
With this information, we find the ratio 



a gg a\ f^\ 2 0.15 / \2 



'66 



3tt 2 tan 4 (3 \m b ) tan 4 (3 \ 100 GeV / 



(3.10) 



Therefore, if the neutralino annihilates primarily through the A , we do not expect the 
gluon final state to be significant. Note that this can be a somewhat artificial example 
since it requires that 2m^ be near vtia for the contribution from the intermediate A to be 
dominant. 

In the limit of large squark masses, the cross section for neutralino annihilation into 
the qqg 3-body final state is suppressed relative to those for annihilation into the light- 
quark and two-gluon final states by a factor of {^^j ■ Although the matrix elements from 
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M x (GeV) 

Fig. 3 Effect of QCD corrections to the cross section for annihilation of higgsinos, photinos, 
and S-inos into light fermions (fr, cc, and bb pairs). 

diagrams (a) and (b) in Fig. 2 are each inversely proportional to m~ in the large squark- 
mass limit, these leading-order contributions cancel and the contribution from diagrams 
(a) and (b) — as well as that from diagram (c) — to the matrix element is proportional to 
m^~ 4 in the large squark-mass limit. Therefore, if the squark masses are large, the cross 
section for annihilation into the 3-body channels is smaller than those for annihilation into 
the gg and qq channels. 

We now illustrate the effect of using running quark masses in the cross section for 
annihilation of neutralinos into fr, cc, and bb pairs. If the neutralino is a pure photino, 
then ajf oc cfrnj-ej-. Quarks have fractional charge, so annihilation occurs primarily 
into fr pairs, and use of the running quark masses decreases the total cross section for 
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annihilation into light fermions only slightly. The same is true if the neutralino is pure 
B-ino. 

On the other hand, if the neutralino is a pure higgsino, the squark-quark-neutralino 
coupling is a Yukawa coupling and is proportional to the quark mass. Therefore, the cross 
section for annihilation into //is proportional to m®, so annihilation occurs primarily 
into bb pairs. If m r b un is the running quark mass, and m ° n - shel1 [ s the on-shell quark 
mass, use of running quark masses decreases the annihilation cross section by a factor of 
(m r q un /m° q n - shel1 ) 6 . This results in a dramatic decrease in the cross section for annihilation 
into quark-antiquark pairs. In practice, unless M2 is extremely large, a neutralino that is 
primarily higgsino will contain some gaugino component, so the actual decrease in the qq 
annihilation cross section will hardly ever be as dramatic as indicated by the pure-higgsino 
case. 

In Fig. 3, we show the effect of QCD corrections to the cross section for annihilation 
into fr, cc, and bb pairs. The ratio of the cross section using running quark masses to that 
using on-shell quark masses is plotted against neutralino masses from 20 to 150 GeV, and 
we illustrate the results for the cases of a pure photino, S-ino, and higgsino. The graph 
suggests that the effect of QCD corrections to the tree-level cross sections is at least about 
10%, and may, in some cases, be as big as 90%. In the next Section we will discuss the 
effect of these results on the energetic-neutrino signals. 

In Fig. 4, we illustrate the importance of the gg and gqq annihilation channels relative 
to the light-fermion annihilation channels for a neutralino of more arbitrary composition. 
In the regions shaded most heavily, the cross section for neutralino annihilation into the gg 
and gqq final states is greater than that for annihilation into light fermions, a gg + a g q q > 
a j j--. In the more lightly shaded regions, the cross sections for the new channels are at least 
0.1 times as big as the cross section for annihilation into light fermions, Ugg+Ugq-q > 0.1 ait. 
In the most lightly shaded regions, the cosmological relic abundance of the neutralino is 
< 0.05, where is the cosmological mass density of neutralinos in units of the 
critical density, and h is the Hubble parameter in units of 100 km/sec/Mpc. In these 
regions of parameter space, the neutralino is viable, but its abundance is too small to 
account for the dark matter in galactic halos. In the empty region below and to the left 
of the solid curve, the mass of the chargino is less than 45 GeV, so this region is not 
experimentally viable. In the empty region in the upper right-hand corner of the graph, 
the neutralino mass is greater than the top-quark mass (which we take to be 150 GeV), 
so our results will have little effect. The short-dash curve indicates the = 80-GeV 
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Fig. 4 Strength of gg and gqq final states relative to the light-fermion final states. The 
darkest region corresponds to a gg + a ggq > crjj, the next lightest to a gg + a gqq > 0.1a ff, 
and the lightest region corresponds to fi x /i 2 < 0.05, as indicated. Also shown are contours 
of = 80,150 GeV and gaugino fraction contours at (0.9, 0.5, 0.1). See the text for 
discussion of other parameters. 

contour, and the diagonal long-dash curves are of gaugino fractions (= + Nq 2 ) of 0.1, 
0.5, and 0.9, as labelled. Therefore, in the bottom right part of the graph, the neutralino 
is predominantly gaugino, and in the upper left part, the neutralino is predominantly 
higgsino. We took tan/? = 2, m A = 500 GeV, m t = 150 GeV, and m q = 200 GeV. We 
also took all squark and slepton masses to be degenerate and assumed no mixing of right 
and left squarks, and we used running quark masses to implement the leading-log QCD 
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Fig. 5 Strength of the gqq final state relative to the gg final state. The darkest region 
corresponds to & qqg /a gg > 1, the next lightest to o qqg jo gg > 0.1, and the lightest to 
Q x h 2 < 0.05, as in Fig. 4. 

corrections. 

Note that the new channels can be quite important for a large range of masses and 
neutralino mixings, and are especially important if the neutralino is primarily gaugino. 
Note that in large regions of parameter space where the neutralino makes a good dark- 
matter candidate, our results are significant. We should also mention that changes to 
Fig. 4 should be small if we change the assumed squark mass. This is because in the 
heavy-squark limit, cr gg /ajf is squark-mass independent. In addition, in the regions where 
we show Vt^h 2 < 0.05, neutralino annihilation occurs primarily into gauge bosons (where 
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the neutralino is heavier than the W-boson and primarily higgsino) or into light quarks 
through gauge and/or Higgs bosons (in the regions where the neutralino mass is roughly 
50 GeV). The cross sections for these processes do not depend on the squark mass. 

In Fig. 5, we illustrate the importance of the 3-body final state relative to the gg final 
state for the same SUSY parameters that were used in Fig. 4. The most heavily shaded 
regions are those where the cross sections for annihilation into gqq are larger than that for 
annihilation into gluons, a g q q > a gg . In the more lightly shaded regions, o gqq > 0.1 a gg . 
Again, the most lightly shaded regions are those where the relic abundance of the neutralino 
is too small to account for the dark matter in galactic halos. Note that the 3-body channel 
seems to be most important only in the regions where the new channels do not have much 
effect (c.f., Fig. 4). If larger squark masses were used, the relative importance of the 3-body 
final state would be even smaller. This suggests that the effect of the 3-body final state is 
generally, although not always, negligible. 

We should point out that, in Fig. 4 and Fig. 5, we illustrated the importance of the 
new annihilation channels under specific assumptions about several SUSY parameters. 
In general, there is a large viable range for several of the parameters, and the relative 
importance of the new channels for a given set of assumptions may be either greater or 
smaller than indicated in the limited regions of parameter space that we have explored. 



4. Energetic Neutrinos From the Sun and the Earth 

If neutralinos populate the galactic halo, then they will be captured in the Sun and 
in the Earth [ IT|j , annihilate therein, and produce high-energy neutrinos that could be 



detected in underground detectors |J 11(| . Neutralinos from the galactic halo are accreted 
onto the Sun and Earth and their numbers therein are depleted by annihilation. Typically, 
the two processes come into equilibrium on a time scale much shorter than the age of the 
solar system, in which case the rate for neutralino annihilation is equal to the capture 
rate divided by 2, C/2. Then, the differential flux of energetic neutrinos of type i (e.g. 
i = 1/^,1?^, etc.) from neutralino annihilation in the Sun or Earth is 

cUf>\ C v-^ D fdN\ 

m)r^^ Br \w» (4 - 1) 

where R is the distance from the detector to the center of the Sun or Earth, the sum on 
F is over all annihilation channels, Bp is the branching ratio for annihilation into channel 
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F, and (dN / dE)pi is the differential energy spectrum of neutralino type i at the detector 
expected from injection of the particles in channel F at the core of the Sun or Earth JT^] . 
Calculation of these spectra can be quite complicated as it involves hadronization of the 
annihilation products; furthermore, if annihilation takes place in the Sun, interaction of 
the annihilation products with the solar medium as well as interactions of the neutrinos 
as they propagate through the Sun must be taken into account. 

In the rest of our discussion, we will focus on the neutrino signal from neutralino 
annihilation in the Earth. If the neutralino is lighter than the top quark, then the neutrino 
signal from the Earth should be greater than or equal to that from the Sun (unless the 
neutralino has only axial interactions with nuclei in which case it is captured in the Sun 
but not in the Earth). Also, calculation of neutrino spectra from annihilation in the Earth 
is much simpler than the calculation of spectra from the Sun. It should be kept in mind, 
however, that our conclusions will also apply to neutrino rates from the Sun. 

The most promising method of detection of energetic neutrinos is observation of up- 
ward moving muons induced by neutrino interactions in the rock below the detector. Given 
the fluxes (d<p/dE)i, the rate for neutrino-induced upward-moving muons may be written 
simply as 

r<ietector = X)A J (Jj^J E 2 dE, (4.2) 

where the sum is over z^ M , which produce muons, and v ^ , which produce antimuons. The 
rate is proportional to some constant Di, times the second moment of the neutrino energy 
spectrum. This is because the cross section for a neutrino to produce a muon is proportional 
to the neutrino energy, and the range of the muon is proportional to its energy, giving an 
overall dependence on the square of the impinging neutrino energy. 

If neutralino annihilation takes place in the Earth, then interactions of the annihila- 
tion products and neutrinos with the Earth can be neglected; thus, the neutrino-energy 
spectrum from a given annihilation channel F can be determined in a straightforward 



way from the results of Monte Carlo calculations |18| . Furthermore, in this case, the en- 
ergy spectra of neutrinos and antineutrinos are the same, so the rate for neutrino-induced 
upward-moving muons from neutralino annihilation in the Earth may be written 



detector 



39 x 10-20 (sft) (B/) 2 E Bf < Wli >* m " 2 yr_1 - (4 ' 3) 



The quantity (Nz 2 ~) F is the second moment of the energy spectrum of neutrinos 
produced from final state F, divided by the neutralino mass squared. Final-state electrons, 
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muons, and u, d, and s quarks will not produce energetic neutrinos. The weak decays of 
r leptons and c and b quarks produce energetic neutrinos, and expressions for (Nz 2 ) 
for these final states have been given by Ritz and Seckel | I8|| . For neutralino masses 
large compared with the light-fermion masses, (Nz 2 )- bb = 0.0195, (Nz 2 ). c = 0.0084, 
and (Nz 2 )- r = 0.0682; thus the fr final state gives the strongest neutrino signal of the 
fermionic final states. If the neutrino annihilates into gauge-boson pairs, then energetic 
neutrinos are produced directly by the decays of the gauge bosons, and (Nz 2 ) w+W _ = 
0.035[l-(m^/4m^ 2 )] and (Nz 2 ) zz = 0.045[1 - (m|/4m^ 2 )] 0. Expressions for (Nz 2 ) 
from final states with Higgs bosons may also be given |§. For annihilation in the Earth, 
the second moment of the neutrino spectrum from a Higgs boson B is roughly (Nz 2 ) B ~ 
(Nz 2 ) , /2, where the sum on / is over all the decay channels of the B boson, and 
is the branching ratio for B decay into final state /. The branching ratios are given in, 
e.g., Ref. |nj. Leading-order QCD corrections to these results may be included by using 
running quark masses instead of the tree-level quark masses. Although we will not need 
it here, (Nz 2 ') for the top-quark final state can be obtained by noting that the top will 
decay predominantly into a b quark and a W boson and using the results for (Nz 2 ) for 
these final states. 

Given Eq. ( fL3|) and the estimates of (Nz 2 ) above, we can see qualitatively the effect 
of including the gg and gqq annihilation channels, as well as the effect of running-quark 
masses on rates for energetic-neutrino events from neutralino annihilation in the Earth. 
The gluons will produce essentially no energetic neutrinos. Energetic neutrinos will come 
from weak decays of heavy quarks, and only a small fraction of the gluon energy goes into 
heavy quarks. Thus, if the cross section for annihilation into gluons is appreciable, the 
branching ratios Bp into the annihilation channels that do produce energetic neutrinos are 
decreased, and the event rate, given by Eq. ([Of), is decreased accordingly. 

The effect of annihilation into gqq (q = 6, c) is similar, though not as severe, as the 
effect of annihilation into gluons. Again, the emitted gluon will not produce any energetic 
neutrinos. The quark and antiquark will produce neutrinos, but their energies will be 
decreased leading to a smaller event rate. In the heavy-squark limit, = m?/m^ 2 — ■> oo, 
the gluon carries away 1/3 the available energy, which suggests that the rate for energetic- 
neutrino events from the 3-body final state containing c or b quarks would be roughly 
4/9 what it would be if the neutralino annihilated into qq. As the ratio is decreased 
towards unity, the gluon carries away a larger fraction of the available energy, thereby 
weakening the neutrino signal. Moreover, more than 50% of the total qqg contribution 
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Fig. 6 Effect of QCD corrections to the tree-level cross sections for \\ —* QQ on energetic- 
neutrino event rates. We illustrate the results for the cases of a pure photino, pure 5-ino, 
and pure higgsino. 

comes from light (u, d, s) quarks. Therefore, the energetic-neutrino flux from the 3-body 
final state is generally small: If ~ 1, then the branching ratio may be significant, but 
the neutrino signal is small; if r<j is large, the branching ratio into the 3-body final state 
is suppressed. We have not done the calculation of the neutrino spectrum from this final 
state more precisely than indicated here. A more precise calculation would simply involve a 
convolution of the neutrino energy spectra from quark-antiquark pairs with the differential 
cross section in Eq. ( |2.17| ). 

Although the new annihilation channels tend to decrease the event rate, leading- 
order QCD corrections to the tree-level cross section for neutralino annihilation into b 
and c quarks have the opposite effect. QCD corrections decrease the cross sections for 
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annihilation into b and c quarks, and therefore increase the branching ratio into f r leptons. 
r leptons provide a stronger signal, so the neutrino event rate is increased. In Fig. 6, we 
illustrate this effect by plotting the ratio of the neutrino event rate obtained using running 
quark masses divided by that obtained using on-shell quark masses. We illustrate the 
results for the cases of a pure photino, pure S-ino, and a pure higgsino. In all three cases, 
we have assumed annihilation occurs only into light fermions. The effect is largest for the 
photino. This is because the cross sections for annihilation of photinos into r leptons is 
comparable to that into quarks. For S-inos, the branching ratio to r leptons is larger, 
and for higgsinos, the branching ratio to the b quark is larger. Recall that the energetic- 
neutrino flux from a given annihilation channel depends on the branching ratio into that 
channel, and not on the total annihilation cross section. 

In Fig. 7, we illustrate the combined effect of including the gg and gqq final states, 
as well as the QCD corrected cross sections for annihilation into qq pairs, on the rate for 
energetic-neutrino events. We consider the ratio of the event rate predicted using our new 
results to that which would be predicted ignoring the new annihilation channels and the 
QCD corrections to the light-fermion annihilation cross section, (Nz 2 )q CD / (Nz 2 ) tree . 
In the most heavily shaded regions (where the neutralino is primarily gaugino), the ratio 
is less than 0.5. In the more lightly shaded regions the ratio is less than 0.9. In the rest 
of the graph, the ratio is greater than 0.9, but it is nowhere greater than 1.1. Again, the 
most lightly shaded regions are those where the cosmological abundance is too small to 
account for the dark matter in galactic halos. The SUSY parameters used here are the 
same as those used in Fig. 4 and in Fig. 5: niA = 500 GeV, m t = 150 GeV, tan/3 = 2, and 
rriq = 200. 

We have included all annihilation channels in this graph. For our numerical work, we 
assumed that the gqq final state produced no energetic neutrinos. This underestimates the 
true flux, but we are confident that, if the calculation were performed more precisely, the 
results in Fig. 7 would not be altered. We believe so because our estimate of the neutrino 
yield from the 3-body final state is small, and because Fig. 4 and Fig. 5 suggest the new 
annihilation channels are important only when the 3-body final state is subdominant. As 
was the case for Fig. 4, Fig. 7 will change little qualitatively if the squark mass is changed. 

The graph illustrates that our results are most important for gauginos, as discussed 
in the previous Section. Also, note that our results are important in large regions of 
parameter space where the neutralino makes a good dark-matter candidate. Again, we 
should point out that we have explored only a restricted region of parameter space. The 
effect of our results on energetic-neutrino fluxes may be larger or smaller depending on the 
specific SUSY parameters used. 
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Fig. 7 Combined effect of new results on energetic-neutrino event rates. The darkest region 
corresponds to (Nz 2 ) QCD / {Nz 2 ) free < 0.5, the next lightest to (Nz 2 ) QCD / (Nz 2 ) tree < 
0.9, and the lightest region corresponds to O x /i 2 < 0.05, as in the previous plots. 

5. Summary and Discussion 

We have calculated the cross sections for annihilation of neutralinos into the two-gluon 
and gluon-quark-antiquark final states in the nonrelativistic limit. We have also calculated 
QCD corrections to the tree-level matrix elements for annihilation of neutralinos into 
quark-antiquark pairs. These new results should have little impact on existing neutralino 
cosmological-abundance calculations, although they will be significant for event rates for 
indirect-detection schemes. 
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Since neutralino annihilation into light quarks and leptons is helicity suppressed in 
the nonrelativistic limit, the rate for annihilation into gluons, although suppressed by a 2 s , 
may be comparable to or greater than that for annihilation into light quarks and leptons. 
Very few energetic neutrinos are produced by hadronization of gluons in the Sun and in 
the Earth, so annihilation into gluons tends to decrease the rate for energetic-neutrino 
events. The new annihilation channels should have a significant effect if the neutralino 
annihilates predominantly into light fermions. If the neutralino is heavier than the top 
quark, or if it is primarily higgsino and heavier than the W boson, then it will annihilate 
predominantly into top quarks or gauge bosons, respectively, and the new channels should 
have little effect. 

QCD corrections decrease the rate for annihilation into light quarks. The flux of 
energetic neutrinos depends on the branching ratios into the various final states, and not 
on the total annihilation cross section. In addition, the neutrino signal from f r pairs is 
stronger than that from light quarks. Therefore, since the branching ratio into f r pairs is 
increased, while the branching ratios into light quarks are decreased, QCD corrections to 
the process XX ~ * QQ increase the energetic-neutrino flux. 

The combined effect of the new annihilation channels and QCD corrections is greater 
than 10% in large regions of parameter space. If the neutralino is primarily gaugino, 
the rate for annihilation into gluons may be greater than the rate for annihilation into 
light fermions, and the flux of energetic neutrinos will be decreased dramatically. In some 
regions of parameter space, the two effects yield a slight increase in the neutrino flux, 
although this increase is no greater than about 10% in any regions of parameter space we 
explored. 

We have not performed a detailed calculation of the flux of energetic neutrinos that 
come from the gqq final state. If the squark is only slightly heavier than the neutralino, 
then the gluon carries away most of the available energy, and the energetic-neutrino yield 
should be small. On the other hand, if the squark mass is large, the branching ratio for 
annihilation into the 3-body final state becomes negligible. 

In conclusion, the gg and gqq final states, as well as QCD corrections to the tree-level 
amplitudes for annihilation into qq pairs, are appreciable for many regions of parameter 
space. These new results should therefore be included in analyses that constrain SUSY 
dark-matter candidates from limits on fluxes of energetic neutrinos from the Sun and 
Earth. 
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The new results will also have implications for cosmic-ray antiproton searches ||. If 
neutralinos populate the galactic halo, they will annihilate and produce low-energy antipro- 
tons. Cosmic-ray antiprotons are produced in standard models of cosmic-ray propagation 
by spallation of cosmic rays in the interstellar medium, and therefore provide a background 
to the signal from WIMPs. However, this background decreases dramatically for cosmic- 
ray antiproton energies less than about 1 GeV. There are many astrophysical uncertainties 
associated with the predicted fluxes of cosmic-ray antiprotons from WIMP annihilation, so 
nonobservation of such cosmic rays cannot be used to eliminate dark-matter candidates; 
on the other hand, under certain reasonable assumptions about the relevant astrophysics 
and particle physics, the predicted flux of low-energy cosmic-ray antiprotons will be large 
enough to be distinguished from background. Observational upper limits to the cosmic- 
ray antiproton flux are currently about an order of magnitude larger than the background 
expected [^DJ. By performing similar cosmic-ray experiments at the South Pole, the sen- 
sitivity of these experiments can be improved by several orders of magnitude pi], so a 
cosmic-ray-antiproton signature for neutralino dark matter in the halo should be observ- 
able, should such a signature exist. 

Our results will have an effect on the cosmic-ray antiproton flux from neutralino 
annihilation in the halo with magnitude similar to that on rates for energetic-neutrino 
events from annihilation in the Sun and Earth; however, the effect is converse to that 
discussed above. QCD corrections decrease the cross sections for annihilation into light 
quarks, so the number of antiprotons produced from hadronization of light quarks will 
be decreased. On the other hand, if annihilation into gluons is appreciable, then the 
antiproton flux will be increased because hadronization of gluons will produce antiprotons 
(although Monte Carlo results will be needed to determine the precise contribution). The 
enhancement will be most significant if the neutralino is primarily gaugino. We should also 
point out that if the new annihilation channels decrease the neutrino rate, then they will 
increase the cosmic-ray antiproton flux. Therefore, cosmic-ray antiproton searches provide 
an excellent complement to energetic-neutrino searches. 
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